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Abstract 

A dilute suspension of impurities in a low density gas is described by the Boltzmann and 
Boltzman-Lorentz kinetic theory. Scaling forms for the species distribution functions allow an 
exact determination of the hydrodynamic fields, without restriction to small thermal gradients or 
Navier-Stokes hydrodynamics. The thermal diffusion factor characterizing sedimentation is identi- 
fied in terms of collision integrals as functions of the mechanical properties of the particles and the 
temperature gradient. An evaluation of the collision integrals using Sonine polynomial approxima- 
tions is discussed. Conditions for segregation both along and opposite the temperature gradient 
are found, in contrast to the Navier-Stokes description for which no segregation occurs. 

PACS numbers: 45.70.Mg, 05.20.Dd 
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I. INTRODUCTION 



Consider a granular mixture of two mechanically different species in a steady state with 
number densities nQ{r) and n{r), respectively. One component is dilute with respect to 
the other, no(r*)/n(r) << 1, such that this component has negligible effect on the host gas. 
Moreover, the latter is at sufficiently low density that the granular Boltzmann kinetic theory 
applies for its intra-species collisions. The dilute component has negligible intra-species 
collisions and its collisions with the host gas are described by the granular Boltzmann- 
Lorentz kinetic theory The objective here is to provide an exact description of segregation 
induced by a temperature gradient in this context. The motivation is the description some 
years ago of an exact solution to the Boltzmann equation for a steady state with constant 
temperature gradient [2], Q]. That analysis is extended here to include the presence of the 
dilute component with a complementary description of the exact solution to the Boltzmann- 
Lorentz equation. Since there is no limitation on the size of the temperature gradient, the 
results given here extend previous results on thermal segregation obtained from the Navier- 



Stokes equation restricted to small gradients [4]. For the dilute conditions considered here, 
and the absence of gravity, no segregation occurs at Navier-Stokes order in contrast to the 
results obtained here. 

The particles of the dilute component will be referred to as the "impurities" . The hydro- 
dynamic fields obtained for the host gas are zero flow velocity, constant temperature gradient 
in the x direction, dT{x) /dx = 6 , and a constant uniform pressure p = n{x)T{x). The impu- 
rities have a temperature profile Tq{x) proportional to the host temperature Tq{x) = •yT^x), 
and a non-trivial density no(x) expressed in terms of the host temperature field. In the 
dilute limit, the concentrations are po{x) — no{x)/n{x) and p{x) = 1 — po{x). They have the 
relationship dp^/dx = —dp/dx so any spatial variation of po{x) implies the opposite vari- 
ation of p{x) and segregation occurs. Here the segregation is induced by the temperature 

gradient, and it is common to introduce a thermal diffusion factor A defined by 

. d\iaT(x) dlnpoix) 

This dimensionless factor depends on the properties of the two components, A = 
A(a, ao, cr/cTo, m/mo, 6^*), where a,ao are the restitution coefficients for the host-host and 
impurity-host collisions, a, do and m, mo are the species diameters and masses, and 9* = 
is the dimensionless temperature gradient, d being the geometrical dimension of the 



system. In principle, A can be positive or negative within this parameter space. The case 
A = imphes no segregation, while A positive (negative) implies the impurities increase 
concentration against (along) the temperature gradient. This is the thermal analogue of the 
Brazil nut and reverse Brazil nut effects for gravitational segregation js-Sj. 

The distribution functions for the two species are of a "normal" form, meaning that their 
dependence o n sp ace and time occurs entirely through the hydrodynamic fields, n{x),T{x), 
and no{x) jo], llol |. Thus, boundary conditions do not occur explicitly but only through the 
determination of these fields. For example, no external driving source is required in the 
kinetic equation for a stationary state, since this is implicit in the time independence of 
the fields. Instead, the stationary form of the fields is determined self-consistently from 
moments of the kinetic equations. This self-consistency also determines the temperature of 
the impurities as being proportioal to the host temperature, To(x) = •yT^x), with 7 7^ 1 
in general. No reference to hydrodynamics is made, although these moment equations are 
equivalent to the balance equations forming the basis for a hydrodynamical description. 

The steady state obtained occurs by establishing a gradient of the heat flux to compensate 
for local energy loss due to coUisional cooling. Thus it is special to granular fluids and links 
the temperature gradient to the degree of inelasticity rather than to boundary conditions. 
This is similar to steady uniform shear flow where the steady state is possible due to a 
balance of viscous heating and collisional cooling, such that the velocity gradient (shear 
rate) is linked to the degree of inelasticity. In both cases, the control needed to assure 
Navier-Stokes hydrodynamics is lost. In the present case, smaller gradients entails smaller 
pressure at constant restitution coefficient, or smaller inelasticity at constant pressure. Such 
non-Newtonian steady states are a characteristic of granular flows and segregation for such 
states can be qualitatively different from that from Navier-Stokes hydrodynamics. This has 
been illustrated recently for thermal segregation under uniform shear flow 

The next section defines the system and its kinetic theory description. In section IIIII 
scaling forms for the distribution functions are introduced and the implications for the 
hydrodynamic fields are obtained. Three constants must be determined self-consistently. 
One of these, the temperature gradient 6 has been obtained in 2|, |3|. Collision integrals for 
the other two are obtained here. The form of the thermal diffusion factor A is given in terms 
of these constants, and the sign of A is discussed based on approximate evaluations of the 
collision integrals given in the Appendices. 
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II. KINETIC THEORY 



Consider a one component gas of smooth, inelastic hard spheres [d = 3) or disks 
{d = 2) with diameter a and mass m at low density. Their distribution for position r and 
velocity v at time t, f{r, v, t), is determined from the Boltzmann equation (without external 
forces) Q 

[9t + v-^^f = J[v\fJ], (2) 
where the collision operator J [v|/, /] is 

J [v|/, /] ^ a'-' I dv, I da Q{a ■ g){a ■ g) [a-^f{T, v', t)f{r, v[, t) - f{r, v, t)f{r, v,, t) 

(3) 

Here g = v — vi is the relative velocity of the colliding pair, B is the Heaviside step function, 
da- is the solid angle element about the direction of the unit vector a, and a is the restitution 
coefficient characterizing the degree of inelasticity (0 < a < 1). The velocities v', v[ denote 
the restituting velocities for the pair v, Vi, 

v' = V -^-^^^—{a ■ g)a, v[ = vi+ '^ (g- ■ g)a. (4) 

Now consider M additional impurity particles in this gas, all the same but mechanically 
different from the fluid particles. For M N, the primary collisions for the impurity 
particles are with the host gas particles, and impurity-impurity collisions and effects of the 
impurities on the gas distribution function / can be neglected. The distribution function for 
the impurities, F{r,VQ,t), is governed by the corresponding Boltzmann-Lorentz equation, 

[d, + vo-^^F = I[vo\FJ], (5) 

where the operator I [vo\F, f] describes changes in F due to binary collisions between the 
impurity and gas particles, 

I [vo\F, f] = cf"-' I dvr J da eia-goi){a-goi) [a^^Fir, v'„ t)f{r, v[, t) - F{r, vo, t)f{r, v,, t) 

(6) 

gQi = Vq — Vi. The restituting velocities v'q, v[ in this case are 

m f 1 + Ofg M mo (l + ^) 
^o = '»^o ^ ^(o--goi)o-, v[ = vi + ^ ^-{a-gQi)a. (7) 

m + niQ m + mo 
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In the above expressions, a = (cr + <7o)/2, and ctq, mo, and ao are the hard sphere diameter, 
mass, and restitution coefficient for the impurity particles, respectively. 

The macroscopic state of this system is described by the fluid number density n{r,t), 
temperature T{r,t), and flow velocity u{r,t), defined in terms of the distribution function 
by 



n{v,t) ^ 




( 1 \ 


ln{r,t)Tir,t) 






n{r,t)u{r,t) j 







f{r,v,t), 



(8) 



with V{r,t) = V — u{r,t). It is convenient to introduce corresponding fields for a macro- 
scopic description of the impurity particles, 

/ ^Jr- \ / 1 \ 



V 



no(r,t) 
ino{r,t)To{r,t) 



1 

V ^0 y 



Fir,Vo), 



(9) 



with Vo(r,t) = vq — uo(r'o,t). Instead of an impurity velocity, the more usual number flux 
notation jo = ^0^0 has been used. 



III. SCALING SOLUTIONS 



In reference a solution to the Boltzmann equation was described for the special case 
of a scaling form in terms of the hydrodynamic variables. 



f{x, v) = n{x) 



m 



-1 d/2 



2T{x) 



m 



nl/2 



2T{x) 



V. 



(10) 



Such a solution, where the space and time dependence of the distribution function occurs 
only through the hydrodynamic flelds, is called "normal" . The definitions of the fields in 
dH]), and the choice of u = give the self-consistency conditions on (c, ) 

dc 



d 
2 

\0/ 



c 



c . 



(11) 



Here, a similar scaling solution for the impurities is sought, 

d/2 



F{x,Vo) = no{x) 



2To{x) 



$(co) 



Co 



2n{x) 



.1/2 



(12) 



The definitions ([9]) then give the conditions on $, 

/ 1 \ (l\ 



1 I ( ma U/2„- , 



■-0 



$(co) 



(13) 



In order for (fT2|) to be "normal" , it should depend only on the hydrodynamic fields for the 
gas and impurities, i.e., on nQ{x), n{x), and T{x). Dimensional analysis then requires that 
To(x) must be proportional to T(x), 



To(x) = 7T(x). 



(14) 



The constant 7 must be determined in course of solving the kinetic equation (as discussed 
below). Further comments on the implications of normal solutions is provided in the last 
section. 

In terms of these scaling solutions and dimensionless velocity variables, the Boltzmann 
and Boltzmann-Lorentz equations become 



I 1 dlnn 



d\nT d 



CQx 



na 



1 dhiriQ 



dx 



■d-l 



<^'(co) 



2na'^~^ dx dc 

1 d\nT d 

2nW^^^ dx dco 



[c0(c)]| = :r[c|0,0], 

■ [co$(co)]l =X[co|$, 



with the dimensionless collision operators 



J[c\(p,' 



dci / dd- Q(d- ■ w){d- ■ w) 



a 



c -0 c 



(15) 

(16) 

(17) 



X[co|$,0] = j dci j dae{a-Wo){a-Wo) X [ao2$(c'Q)0(c;) -<l>(co)0(ci; 
The relative velocities w and Wq are now 



W = C — Ci, Wq = Cq 



mo 



1/2 



c. 



(19) 



The expressions of the dimensionless restituting velocities in Eq. ( !T8|l are given in Eq. flASp . 
Since the right sides of Eqs. (fT5|) and ( IT6|) are independent of x, the left sides must be as 
well. This will be true if the hydrodynamic fields n{x), no{x), and T{x) satisfy the equations 



1 dlnn 



<^-^ dx 



A, 



dlnT 



na 



d-l 



dx 



B 



d\\i no(x) 



dx 



(20) 
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where A, B, and C are constants. The constants A and B are determined by taking moments 
of the Boltzmann equation f|T5l) . Namely, multiphcation of the equation by 1, c^, and c^, 
and integration over c yields 

[a + ^^J dc c,(t) (c) = 0, {A + B) j dc cl<P (c) = 0, (21) 

(a + ^) / dc ^c^<P (c) = j dc c^J [c|0, 0] . (22) 

The zeroes on the right sides of (1211) result from conservation of particle number and momen- 
tum by the collision operator. The first equation of ( 12T|) is satisfied because of conditions 
f|TT]) required on 0(c), while the second equation gives A = —B. Finally, Eq. fl22|) determines 
B, 

_ 2jdcc'J[c\4>,<l>] 

Jdcc^c^ic) ■ ^'^^ 

The exact hydrodynamic fields for the gas are now given exactly by n = and 

where p = n{x)T{x) is the uniform pressure, and 6 = Bpa'^~^ is the constant temperature 
gradient. 

A similar analysis applies for the impurity constants C and 7. Taking moments of the 
Boltzmann-Lorentz equation f|T6|) with respect to l,Cox, and Cg gives 

(c + h, = 0, (25) 

(T"'-l JdCocl,.<^{co) 

3a''-' /^CocgX[co|$,0] 
'^ + 2a'^~i'' /rfcoc2co.<l>(co) " ^^'^ 
The right hand sides of Eqs. fl2Bl) and flTTl) depend on 7 explicitly through Wq (see f lT^ ) 

and implicitly on both 7 and C through $. Since i? is known independently from Eq. (l23ll . 
the two unknowns 7 and C are determined by Eqs. ( 126|) and (!27|) . Equation (!25|) has two 
solutions, jox = and C = -a'^-'B /2a''-\ The latter gives an additional equation for 
7 and C and the problem is overdetermined. Probably, this choice is not consistent with 
the assumption ( |T2l) . Here, it is assumed that the boundary conditions enforce the choice 
3 Ox = 0. 



7 



In summary, the description of the gas and impurities is completely specified by the 
kinetic equations for (c) and $ (cq), 

1 d 



-BcA<P{c) + -—-[c<P{c)] 



Cox \ C$ (co) - 



J[c\(l),(l)], 
[co<I>(Co)]| =X[co|$,0] 



(28) 



(29) 
, (I2H1), and 



and the constants B, C, and 7 are determined self-consistently from Eqs. 
( 1271) . The corresponding collision integrals are further simplified in Appendix |Al The 
hydrodynamic fields have the simple spatial forms 

P 



T{x) = T(0) + ex, n{x) 



To(x) = To(0) + 76'x, no{x) = no(0) 



T(0) + 9x ' 

1 + 



ex 



(30) 
(31) 



IV. SEGREGATION 



The segregation of impurity particles relative to the host gas is described by the inhomo- 
geneity of the composition po{x) ~ nQ{x)/n{x), which follows from fl3Up and 0311) 

-1 l+W'^-'^C/a'^-'^B 



Po [X) 



nojO) 
n(0) 



1 



e 



r(o) 



X 



The thermal diffusion factor of ([T]) is therefore 



A = - 1 + 



/ dc(?Cx<\) (c) 



(32) 



(33) 



V cy'-^B} 2a'^-i /rfcoc2,.$(co) \ dac^J y\, , ^ ^ 
If the impurities are mechanically equivalent to the host particles, then X[co = 
J7'[c|0, 0], $ = 0, and A = 0, since the first integral in the numerator of (1331) vanishes 
by conservation of momentum. 

The corresponding result for the thermal diffusion factor obtained from the Navier-Stokes 
order Chapman- Enskog solutions to the Boltzmann and Boltzmann-Lorentz equations gives 
A = for all values of the parameters a, ao, c/ctq, m/mo, 6'/pcr'^~^. If the Navier-Stokes 
calculation is extended to include effects of gravity the condition becomes {4] 



, dx , 



NS 



/To 

mg\ — 

I m 



(34) 



Thus thermal segregation can occur, facilitated by gravity, and depends on the sign of 
{Tq/T — mo/m) and the direction of dT/dx relative to the gravitational force. This is in 
sharp contrast to the results obtained in the next section. 
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V. APPROXIMATE DETERMINATION OF Tq/T AND A 



To determine the coefficients B, C, and 'y = Tq/T, the distribution functions and $ are 
represented as truncated Sonine polynomial expansions 



~ TT ' e 



aoi (c^ - dcl^ + (^^-—^boi + ^ho) 



-boic^Cx - (610 - boi) cl 



(35) 



-BqicIcqx — (-B10 — Bqi) C\ 



Ox 



(36) 



The method for determining the coefficients in these expansions is described in [2| and 
summarized for the case here in Appendix [Bl The numerical solutions for the case of a 
two-dimensional system [d = 2) with m = mo and a = ao = are shown in Figs. [71|9] as a 
function of a for several values of Oq- An important general feature is that all coefficients 
in fl35|) and f l36|) vanish as a — )■ 1. Thus the non-uniform steady state described here exists 
only as a consequence of the inelasticity of the host gas. Further comment on this is given 
in the last section below. 

In the following, attention is restricted to cr = ao = a" and d = 2 for several values of 
mo/m, a, and Oq. It is well-established that different species of granular mixtures have 
different partial temperatures, eyen in their homogeneous cooling state (i.e., equipartition of 



energy does not occur) 



14| . Figures [T] and [2] show the behavior of Tq/T for mo/m = 1 



and 2, respectively, as a function of a for several values of Oq. The common feature is 
increasing Tq/T with decreasing a, increasing a^, and decreasing mg/m. Figure E] shows a 
broader range of mo/m. Even for the relatively weak dissipation values of this figure, it is 
clear that the largest values of Tq/T occur for small mass ratio, maximum host dissipation, 
and weakest impurity dissipation. 

The existence of segregation for the same weak dissipation values of Fig. |3]is demonstrated 
in figure m The thermal diffusion factor A is positive for m^/m > 1. This means that the 
impurity concentration is higher at the colder part of the host fluid. This is similar to 
the host fluid density which behaves as n = p/T with constant p. For smaller mass ratio, 
segregation goes in the opposite direction with the impurity concentration highest in the 
hotter part of the host fluid. This effect is enhanced at stronger host fluid dissipation 




0.5 I ' ' ^ 

0.4 0.6 0.8 

a 



FIG. 1: Temperature of the impurity To divided by the temperature of the hot gas T as a function 
of the coefficient of normal restitution of the gas particles a, for several values of the restitution 
coefficient for collisions between the gas particles and the impurities, ao, as indicated in the insert. 
In all cases, d = 2, mg = m, and a = aQ = W. 




0.8 - 



0.4 I ' ' ' ' 1 

0.2 0.4 0.6 0.8 

a 

FIG. 2: The same as in Fig. [H but now mo/m = 2. 

and weaker impurity dissipation, as illustrated in figures E] and for mo/m = 1 and 2, 
respectively. It is interesting to note that for mo/m = 1 the border between the two types 
of segregation, A = 0, occurs for a = ao. Referring to figure HI these values also correspond 
to Tq/T = 1. Similarly, for mo/m = 2 comparing Figs. [2] and |6l it is seen that A = for 
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1.2 



1 

0.8 



0.6 
0.4 

2 4 6 8 10 

rrio/m 

FIG. 3: Temperature of the impurity To divided by the temperature of the hot gas T as a function 
of the mass ratio mo/m for several values of the coefficients of normal restitution a and ao, as 
indicated in the insert. In all cases, it is d = 2 and a = aQ. 

To/T = 2. This limited data suggest the possibihty that the segregation criterion A = 
occurs for mo/m = Tq/T. Surprisingly, this is the same as the Navier-Stokes criterion in 
the presence of gravity, fl34|) . Further analysis of this potential relationship across a larger 
data set is required. For larger mQ/m it is found that Tq/T < 1, and only the segregation 
for A > occurs. 




VI. DISCUSSION 

The description of a low density granular gas with a dilute concentration of impurities 
has been given in terms of solutions to the coupled Boltzmann and Boltzmann-Lorentz ki- 
netic equations. These are normal solutions whose space and time dependence are entirely 
specified in terms of the hydrodynamic fields n, no, and T. The special case of a steady state 
in which the host gas has a constant temperature gradient and constant pressure, described 



earlier in refs. 3 and fl, has been generalized to include a corresponding steady state 
of the impurities. In this way the thermal segregation factor is identified in terms of the 
constants of the hydrodynamic fields, without the limiting approximations of small spatial 
gradients. The self-consistent kinetic equations and (12^ determining these constants 
was solved using a low order Sonine polynomial approximation for the velocity dependence 



11 



1.5 
1 

< 0.5 


-0.5 

2 4 6 8 10 

rrio/m 

FIG. 4: Dimensionless thermal diffusion factor A as a function of the mass ratio mo/m for the 
same system as in Fig. [3l 

0.5 
0.3 



< 0.1 
-0.1 
-0.3 

0.4 0.6 0.8 

a 

FIG. 5: Dimensionless thermal diffusion factor A as a function of the coefficient of normal restitution 
of the gas particles a, for several values of the restitution coefficient for collisions between the gas 
particles and the impurities, oq, as indicated in the insert. In all cases, d = 2, rriQ = m, and 
a = ao = W. 

of the host and impurity distributions. The resulting thermal diffusion factor was found to 
identify conditions for both segregation along and against the temperature gradient. Such 
normal solutions are typically constructed by the Chapman-Enskog method whose practical 
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FIG. 6: The same as in Fig. [U but now rriQ/m = 2. 



application typically entails limitations to small spatial gradients, e.g. Navier-Stokes order. 
Application of Navier-Stokes hydrodynamics obtained in this way, and specialized to the 
steady state with constant temperature gradient and constant pressure, leads to the predic- 
tion of no segregation. The effects described here therefore are due to contributions from 
the Chapman- Enskog method beyond the small gradient approximation. In fact, there are 
no limitations on the temperature gradient in the present analysis. 

There are two important clarifications to note. First, the validity of a normal solution 
both for granular and molecular gases is limited to domains away from the initial preparation 
time and confining boundaries. For the steady state considered here, this means that there 
is typically a boundary layer across which the normal solution does not apply. Additional 
information is then required to connect the physically specified values of the fields or their 
gradients at the boundary with those values associated with the normal solution. These are 
the familiar " slip" boundary conditions. The existence of the normal solution described here 
for a system with finite confinement and associated boundary layer has been demonstrated 
by molecular dynamics simulation in refs. Q and fl. Typically, the size of the bulk interior 
relative to the boundary layer decreases as the temperature gradient is increased. Investi- 
gation of this problem for a molecular gas has demonstrated that the bulk normal solution 
domain still exists beyond the Navier-Stokes limit jl5 |. 

A second clarification is the special nature of the steady state described here as being 
unique to a granular gas. The analysis of shows that it results from the balance of the 
heat flux gradient and the cooling rate due to inelastic collisions. In the absence of the 
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latter there is no steady state solution of the type considered here. In contrast to normal 
fluids, the gradients of such steady states are controlled by internal processes rather than 
boundary sources. External control of the gradients is therefore lost. In the present case 
the magnitude of the dimensionless temperature gradient Qj-po^^^ = B (a) monotonically 
decreases to zero as a — )■ 1, vanishing in the elastic limit. Consequently, for example, it is 
not possible for the Navier-Stokes to apply here for strong dissipation. 
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Appendix A: Reduction of collision integrals 



The Boltzmann collision integral appearing on the right side of Eq. f l2^ is simplified 
further in 2|, with the result 

(1 - a^) Tji'^-^y^ J dcjdci \c -cifcf) (c) (ci) 



B 



(Al) 



2r(^^) Jdcc^c^cPic) 

The Boltzmann-Lorentz collision integrals can be simplified in a similar way. Consider first 
the collision integral appearing in Eq. 



dco cq^I [co|$, 4>] = J dci J dco cq^ J da Q{a-Wo){a--Wo) [a^ $ (cq) {c[) - $ (cq) (p (q) 

(A2) 

where Wq is defined in Eq. ( [T9|) and the dimensionless restituting velocities following from 
Eq. ([7]) are 



m 



Cq = Co 



m + mo 
It is easily verified that 



1 + c^o ""^j ■ "^o)^, c = c 



—- ( 1 + "o ) ( ) ^ ■ "^0)0-. 

m + mo ^ ^ mo 



c/c'^Cq = tto dcidco, 5- ■ go = -aoa ■ gg. 



Also, Eqs. (IA3p can be inverted to get the collision rule in dimensionless units, 

m (1 + ao) 



Co = Co 



m + mo 
14 



(A3) 



(A4) 



(A5) 



m + rriQ \ ttlq 

Returning to Eq. (1A2I) . change variables in the first term of the brackets on the right hand 
side to integrate over the restituting velocities. Using the above relations, the equation 
becomes 

j dcoCoxX[co|$,0] = J dci J rfco $ (co) (ci) J da Q{a- ■ wo){a- ■ Wo) {c*Q^ - cox) 

m (1 + ao) J J ^ ^^^^ ^ ^^^^ wqWox J da- 6(5- ■ Wo){a- ■ wof 



m + niQ 
m(l + ao) 7r('^"i)/2 



m + mo Y 

Finally ( 126|) becomes 



d+3 
2 



J dci J dco^{co)(l){ci)woWox. (A7) 



^ ^ (^'^ m (1 + Qq) TT^'^-^y^ J dCi J dCp $ (Cq) (f) (Ci) WqWqx 

m + mo T (^) / c/cq cg,<l> (cq) 

The analysis of Eq. (!27|) is similar with the result 



m (1 + tto) Tr^'^-i)/^ / dci J dco $ (cq) (ci) 




m + mo r 


/rfCoC§Coa;$(Co) 



Appendix B: Solutions to kinetic equations 

The solution to the kinetic equation for (c) and the self-consistent determination of B 
is a problem that is independent of the impurities and can be carried out first. The method 
is described in First, B is given its representation as a collision integral using Eq. 
so the kinetic equation fl25]) becomes 



Next (c) is approximated by a truncated Sonine polynomial expansion 
0(c) ~ n-'^/^e-'' 



1 - aoi (c^ - dcl^ + ( ^^-^ ^01 + 1^10 ] Cx 



2 

-&oic'c,- (&io-&oi)4j • (B2) 

This form assures the conditions given in Eq. ( ITTj) . The coefficients aoi,6oi, and bio are 
then obtained from three equations following by taking velocity moments in ( IBip . Namely, 
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the equation is multiplied by c^, c^, and c^c^, respectively, and afterwards integrated over c. 
With these coefficients determined, B is calculated from Eq. f l23|) . 

To determine $(co), C, and 7, a similar procedure is followed. First, express C as a 
collision integral from Eqs. f l2B]) and (I27p . 

3 / C/Co Co:r2^ [co I 0] / rfco cgX [co I $, 0] 



(B3) 



2 /c/coc§^$ (co) /c/cocgcox$(co) 

and use this in the kinetic equation (!29|) . Next, express $ (cq) as a truncated Sonine poly- 
nomial expansion 



$(co) 



v4m ( Cn - (iCr 



01 [^0 - ^^Ox ) + 
-BqicIcqx — (Bio — Boi) Coa; 



'd 



(B4) 



which satisfies the conditions (fT3l) with jo = 0. The coefficients, Aqi, Bqi, and Sio are deter- 
mined from three equations obtained by taking moments of fl29|) with respect to Cq^., c^^, and 
coxCq- However, these equations also depend on 7, so they are supplemented by an additional 
equation relating the above coefficients to 7. It is obtained from a new combination of Eqs. 
(EED and (EZD 



B 



2W 



d-l 



J dcoclI[co\^,(t)] J dcoCo^I[co\^,> 



J dcoclcox^ico) Jdcocl^^{co) ' ^^^^ 
Since and B are known at this point, this gives four independent equations for the coeffi- 
cients Aqi, Bqi, BiQ, and 7. With these determined, C is calculated from Eq. (IBSp . 

In practice, the above procedure leads to highly nonlinear equations for the coefficients. 
In the numerical results to be presented in the following, only terms up to second degree in 
the coefficients have been kept sj. As an example, in Figs. UM\ the parameters obtained 
for a two-dimensional system {d = 2) with m = mo and a = ao = a are plotted as a 
function of a for several values of Oq- For small values of a, the numerical solutions for the 
B parameters constructed as described above seem to disappear. 
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FIG. 7: The dimensionless parameters aoi and Aqi as a function of the coefficient of normal 
restitution of the host gas particles a, for several values of the coefficient of restitution for the 
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